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ABSTRACT 


The  classical  solution  for  linear  estimation 
of  a  signal  In  additive  noise  Is  given  by  the  well- 
known  Wiener  filter.  In  certain  circumstances, 
however.  It  Is  advantageous  to  have  a  filter  which 
la  piecewise  constant  and  which  provides  good  per¬ 
formance  over  an  entire  class  of  signal  and  noise 
power  spectral  densities.  First,  a  procedure  Is 
presented  for  determining  the  optimum  piecewise 
constant  filter.  Then,  for  a  particular  class  of 
power  spectral  densities, termed  the  p-polnt  spect¬ 
ral  class,  Che  robust  filter  Is  determined. 

I.  INTRODUCTION 

In  the  linear  estimation  of  a  random  signal  In 
additive  noise,  an  estimate, s (t) ,  of  the  signal, 
s(t),  is  obtained  by  passing  the  received  signal, 
*(t)  “  8(0  +  n(t) ,  through  a  linear  filter  with 
Impulse  response,  h(c,T).  Assuming  the  signal  and 
noise  processes  are  real,  zero-mean,  uncorrelated 
and  wlde-sense  stationary  with  power  spectral  den¬ 
sities  (PSD's),  S(w)  and  N(w),  respectively,  the 
Wiener  filter  solution  Is  obtained  by  choosing 
h^CO,  Che  Impulse  response  of  Che  optimum  filter, 

to  minimize  the  mean-square-error  (MSE)  between  the 
signal  and  Its  estimate.  The  Fourier  transform, 
H^(w),  of  h^(t)  is  Che  solution  Co 

e(S.N;H^)  -  “g“e(S,N;H)  (1) 

where  the  MSE,  e(S,N:H},  Is  given  by 

elS,N;H)  -  ij*iS(w)Il-H(w)I^4«(w)|H(w)[2]dw  (2) 


This  la  for  the  case  where  H^Cw)  Is  not  constrained 
to  be  causal. 

The  solution  to  (1)  Is  found  to  be 


II  a  -  ^  >  t- 

“o'  >  S(w)-Hl(w) 


the  well-known  Wiener  filter.  Knowing  the  signal 
and  noise  PSD's,  one  can  thus  obtain  the  optimum 
Wiener  filter. 

Section  II  details  a  procedure  for  obtaining 
a  piecewise  constant  'filter  which  minimizes  the 
HSE.  By  deriving  the  necessary  conditions  for  min¬ 
imization,  the  optimum  piecewise  constant  (OPC) 
filter  is  shown  to  be  distinct  from  the  optimum 
quantization  (OQ)  of  H^(v).  The  OPC  filter  pro¬ 
vides  a  simplified  filter  structure  with  only  a 
small  degradation  In  performance. 


In  Section  III,  we  show  that  a  piecewise  const¬ 
ant  filter  also  possesses  the  property  of  robustness 
for  a  particular  class  of  PSD's  -  the  .p-polnt  spec¬ 
tral  class.  A  proof  of  the  robustness  of  Che  piece- 
wise  constant  filter  Is  given  along  with  a  justifi¬ 
cation  based  on  distance  measure  concepts.  Several 
authors  have  previously  considered  aspects  of  robust 
filtering.  Kassam  and  Llm  [1]  and  Kassam,  Llm  and 
Clmlnl  [2]  considered  robust  Wiener  and  matched  fil¬ 
ters  using  a  band-model  for  the  PSD's.  Poor  (3,  4] 
Chen  considered  these  filters  In  a  more  general 
framework. 

The  results  obtained  here  are  useful  in  situa¬ 
tions  where  filtering  is  implemented  digitally  (e.g. 
using  FFT's).  In  image  processing  applications,  for 
example,  Wiener  filters  have  often  been  used  [S]. 
Also,  piecewise  constant  filters  may  be  particularly 
useful  in  situations  where  Che  filter  parameters 
must  be  periodically  adjusted.  This  Is  due  to  the 
relatively  simple  forms  of  these  filters. 

II.  OPTIMUM  PIECEWISE  CONSTANT  FILTER 


As  stated  in  the  Introduction,  In  certain  sit¬ 
uations  it  may  be  advantageous  to  use  a  piecewise 
constant  filter.  For  example,  this  is  particularly 
helpful  when  the  filter  parameters  must  be  periodi¬ 
cally  adjusted.  In  this  section,  we  provide  a  pro¬ 
cedure  for  obtaining  the  optimum  piecewise  constant 
filter. 

Considering  symmetric  PSD's  such  that  S(w>'S(-w) 
and  N(w)-N(-w)  with  total  signal  and  noise  powers, 

Og  and  o^,  given- by 


ij 

Ps(w)dw 

(4a) 

pN(w)dw 

(4b) 

we  assume  a  fixed-form  filter  that  Is  piecewise 
constant,  l.e., 

H(w)  ■  Cj,  w  <  1-1, 2,... m  (5) 

where  the  C^  are  real,  non-negative  constants  and 

{a,)“  Increasing  with  a  •  0  and  a_  •  n  where  Q 
J  J»o  “  o  m 

Is  the  known  bandwidth  of  the  signal.  The  result¬ 
ing  MSE  becomes 

e(S,N;H)  -  i  |  ^  [S(w)  (1-Cj)^-Hl(w)cjldw  (6) 

The  breakpoints,  Sj,  and  the  levels,  Cj,  then  must 
be  chosen  to  minimize  (6) . 
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To  alnlnl7,e  the  MSE  we  differentiate  It  with 
tospcct  to  and  Cj.  The  partial  derivative  with 

Toapact  to  Cj  gives 


I-2(l-Cj)S(w)+2CjN(w) Idw. 


Interval  [aj^_j^,  a^) ,  ijlJ),  If  there  exists  a  subdi¬ 
vision  with  different  levels  which  results  In  a 
lower  HSE,  then  the  original  filter  Is  not  optimal. 
The  MSE  associated  with  the  frequency  Interval  [a,b) 
with  filter  level  Is  given  as 

a(S.M;H),^  -  ^lh\^  (W) 


Setting  this  equal  to  zero  and  solving  for  C. ,  we 
obtain  ^ 


i  S(w)dw 

•.a  a 


1 

(S(w)-HI(w))dw 
J-1 


O) 


lAtch  la  simply  the  ratio  of  the  signal  power  In 
the  Interval  [aj_j^,a^I  to  the  corresponding  sig¬ 
nal  plus  noise  power.  The  partial  derivative  with 
zaapect  to  a^  gives 

“  ^<l-Cj)'s(aj)-a-Cj^l)'s(aj)]  + 


Jccj  »(aj)-cj^^I«(aj)]. 

Setting  this  equal  to  zero,  we  obtain  the  equation 


8C.j)  I  (2-Cj-Cj^^)  (Cj^^-C^l-H(aj)  [ 


1  1  J  S  n-l  (8) 

(where  we  are  also  assuming  that  S(w)  and  N(w)  are 
continuous  at  w  -  a^,  1  a  J  in).  Equations  (7) 

end  (8)  are  a  set  of  (2ii»-l)  equations  which  must  be 
satisfied  for  a  piecewise  constant  filter  to  pos¬ 
sibly  be  the  OFC  filter.  Vhen  considering  a  filter 
for  which  for  some  J")^,  then  (8)  Is  always 

satisfied  for  for  other  cases,  (8)  is  equl- 

wslant  to 


■o<‘j>  ■  21 

(provided  we  allow  H  (a.)  to  have  any  value  when 
S(Sj)-ll(a^)-0).  ®  J 

As  stated  above,  when  for  some  J,  (8) 

Is  alwmya  satisfied.  However,  as  the  following  dls- 
•CHSslon  shows,  given  the  number  of  allowable  break¬ 
points,  this  situation  Is  not  optimal  unless  H^(w) 

Is  Itself  a  piecewise  constant  filter.  Consider  a 
piecewise  constant  filter  with  Cj~Cj^^(for  some  J, 

and  such  that  the  necessary  conditions  above  are 
satisfied).  For  any  constant  Interval  (a,b)  with 
lavsl  (l.e.,  either  other 


where 


and 


Sab 

I 

2  .  1  I 

».h  'J 


S(w)dw 


H(w)dw 


With  defined  as  Ir.  (7),  the  partial  HSE  of  CLO) 
becomes 


2  2 
®S  ®ll 

*(S*«!S>ab  “ 

\b 


(11) 


Mow,  using  a  different  filter,  H*(w),  obtained  by 
dividing  the  Interval  (a,b)  into  [a,t)  and  [t,b) 
with  levels  determined  from  (7) ,  with 


»“(t) 


•! 


S(w)dw 


H(w)dw 


(thus  erg(b-t)“Og 

ab  ab 

ve  obtain  the  partial  MSE 


Og(t)oJ(t)  a^(b-t)a^(b-t) 

“^^**''“‘>ab-  izdrz.  *  J::-  2  - 


(12) 


ag(t)-w‘(t)  Og(b-t)-w‘(b-t) 


Ey  writing  out  the  MSE  expressions  and  performing 
some  algebraic  manipulations  we  obtain  the  condi¬ 
tion  for  lowering  the  HSE.  If  there  exists  a  point 
te(a,b)  such  that 


Og(t)  ^  ag(b-t) 


Oj,(t)  Ojj(b-t) 


(13) 


then  the  filter  with  is  not  optimum.  On 

the  other  hand.  If 


for  all  tE(a,b) 


(J*(t)  o2(b-t) 


I'lia  niiits 

DDC  Buff  SecUoQ 

UNANNOUNCCO 

JUS' nON 


Sectkn  |f| 


then 


□  □ 


for  all  te(a,b) 


(15b) 


S<t) 

H(t) 


ab 


"ab 


or,  equivalently,  H^(w)  la  constant  over  the  In¬ 
terval  (a,b).  Thus,  If  (13)  Is  not  satisfied  for 
all  the  Intervals,  H  (w)  Is  a  piecewise  constant 
filter. 

Tor  the  special  case  where  for  all  J, 

the  {liter  Is  a  constant  for  all  w  and  Is  given  as 


H(w) 


2  2 


This  filter  structure  has  previously  been  used  (In 
two-dlnenslonal  form)  In  Image  processing  applica¬ 
tions  [5]  when  very  little  Is  Icnown  about  the  sig¬ 
nal  and  noise  processes. 


-  Sj,  we  oust  simultaneously  solve  (7)  and  (8)  for 

all  j.  In  all  but  extremely  simple  cases  these  e- 
quatlons  are  not  directly  solvable  but  require  an 
Iterative  technique  for  solution.  Max  [6]  devel¬ 
oped  an  algorithm  for  Iteratively  solving  equa¬ 
tions  such  as  these.  For  the  examples  which  we 
will  consider  a  simpler  searching  algorithm  for 
finding  the  minimum  of  a  function  of  several  vari¬ 
ables  was  applied  to  the  general  expression  for 
the  MSE  (6) . 

Equations  (7)  and  (8)  are  only  necessary  con¬ 
ditions.  In  general,  the  sufficient  conditions 
are  difficult  to  obtain.  The  best  answer  is  that 
if  all  the  second  partial  derivatives  of  the  MSE 
with  respect  to  the  levels  and  the  breakpoints  ex¬ 
ist,  then  a  critical  point  determined  by  the  nec¬ 
essary  conditions  Is  a  minimum  if  the  matrix  whose 
ith  row  and  jth  column  element  Is 


I 

critical 

point 


Xj-  Cj  or  Sj,  J-l,2,...m 


is  positive  definite  [6]. 

At  this  point  It  Is  Interesting  to  consider 
the  optimum  quantization  (OQ)  of  H  (w).  The  OQ  of 
R  (w)  Is  obtained  by  minimizing  ° 


f  (3  2 

*  -  I  J  {H„(w)-C^)*dw  (M) 

J-1  Aj.1 

Minimizing  as  before,  we  obtain  the  necessary  con¬ 
ditions 


f  ^  H^<w)dw 

svarage  value  of  H_(w)  on  the 

g  .  -  o 

3  dj  ■  •j.i  interval  taj_j,  a^l 


which  SMist  be  solved  simultaneously  to  obtain  the 
OQ  of  H^(w).  Obviously,  the  OPC  filter  and  the  OQ 

of  H  (w)  are  not  the  same,  and.  In  a  linear  filter- 
o 

Ing  sense,  the  OQ  filter  Is  not  optimum.  That  is. 

It  Is  not  the  best  way  to  obtain  a  piecewise  con¬ 
stant  filter  given  the  signal  and  noise  PSD's.  If 
H^(w)  Is  already  a  piecewise  constant  filter,  the 

OPC  and  OQ  filters  will  coincide. 

As  an  Illustration  of  the  design  of  an  OPC  fil¬ 
ter,  we  will  consider  the  specific  class  of  signal 
and  noise  PSD's  possessing  the  general  Butterworth 
form  (tbM.)  i.e.. 


8(w> - V-2P  .  P  i  1 


(16a) 


and 


H(w) 


B 


l+(f->^ 


,  n  t.  1 


(16b) 


iriiere 


sin 


2P 


and 


Blo¬ 


using  these  as  our  nominal  densities  we  determined 
the  OPC  filter  parameters  and  calculated  the  resul¬ 
ting  HSE's.  Table  1  gives  the  results  for  several 
values  of  p,  n,  and  a-.  In  all  cases,  we  assumed 
2  2'  ° 

Og"Ojj“l  and  ajj“l.  The  MSE  using  the  optimum  filter 


has  also  been  tabulated  for  comparison.  These  re¬ 
sults  show  the  small  degradation  In  performance  when 
using  the  OPC  filter  as  opposed  to  the  optimum  fil¬ 
ter,  Ho(w).  In  each  case,  the  optimum  filter  is 

very  nicely  approximated  by  the  OPC  filter.  That 
is,  the  OPC  filter  provides  the  benefits  obtained 
from  quantizing  H^(w)  but  with  a  smaller  MSE  than 

If  we  were  to  use  the  OQ  filter.  One  should  note 
that  for  this  particular  set  of  PSD's  very  little 
improvement  In  performance  Is  obtained  by  Increas¬ 
ing  m  (I.e.,  finer  "quantization"). 

As  an  example,  we  will  consider  the  case  p-3, 
nwl,  and  a  >1,  i.e.. 


S(w) 


6.i4 


aad 


R(w)  - 

1+w^ 

The  optimum  filter  (  3)  la  given  as 


<15a) 


V 


I 


H  (w) - i - g - 

1+ _ iis! _ 

3  aln|<l-tv^) 

with  the  corresponding  MSE,  e(S,N;H  )-0.393A. 

o 

OFC  filter  (w"!)  Is  determined  to  be 


The 


®OPC<’'> 

(■-2) 


0.615. |w|<  1.392 
0.080, |W|^  1.392 


with  the  corresponding  MSE,  e(S,N:H._.  )-0.4030 

OPC^j 

which  Is  a  degradation  In  performance  of  only  2.4C. 
When  an  additional  breakpoint  Is  allowed  (l.e. , 
■^3),  the  OPC  filter  is  found  to  be 


"opc 

(■-3) 


(w) 


0.627, |w|<  1.1844 
0.352,  1.1844  <Iw|<  1.6807 
0.039. IwU  1.6807 


with  the  corresponding  MSE,  e(S,M;H 


bPC. 


)-0.3972, 


'n«3 

a  degradation  from  using  H^(w)  of  less  chan  IZ.  As 

shown  In  Figure  1,  the  OPC  filter  nicely  approxi¬ 
mates  the  optimum  filter,  H^(w),  at  the  cost  of 

only  a  small  degradation  In  performance. 


III.  ROBUSTNESS  AMD  THE  p-POIMT  SPECTRAL  CLASS 


In  order  to  Implement  the  Wiener  filter,  the 
signal  and  noise  PSD's  must  be  known  exactly.  This 
is  usually  an  unreasonable  assumption  and,  conse¬ 
quently,  small  deviations  from  these  assumed,  or 
nominal,  PSD's  may  result  in  large  degradations  In 
performance.  In  situations  such  as  these  we  would 
like  to  design  filters  which  are  less  sensitive  to 
deviations  from  the  assumed  PSD's. 

Essentially,  we  are  Interested  in  obtaining  a 
filter,  Bj^(w) ,  which  will  give  a  non-trlvlal  upper 

bound  for  the  MSE  over  Che  classes  of  signal  and 
noise  PSD's,  ,i  and  N.  Such  a  filter  Is  robust  In 
Che  sense  that  a  certain  level  of  estimation  per¬ 
formance  Is  guaranteed  for  the  entire  classes,,!  and 
K.  Mathematically,  the  robust  Wiener  f liter ,  Hj^(w), 

Is  defined  Co  be  the  minlmax  filter  such  that 


min  max  max 

SeX  e(S,N;H)  -  SeX  e(S,N;lL)  <17) 

H  NeN  NcK 

If  Is  also  the  optimum  filter  for  a  least  fa- 

worsble  pair  of  PSD's,  S^(w)eX  and  Nj^(w)eH,  Chenwe 

have 

a(S.N:R^)  1  e(Sg,Nj^;Hj^)  s  e(Sg,N^;H)  (18) 

for  any  pair  (S,N)elxN  and  any  linear  filter  R(w). 

The  particular  classes  of  PSD's  which  we  will 
consider  are  defined  as 


>17/*^  S(w)dw  - 

•j-1 


P8*8' 


J-l,2.3....a)  (19e) 


M  -{N(-)l7  N(w)dw  -  p  o^,  j-l,2,3,...m}  (19b) 

•3.1  i 

where  a  ^^0,  a_^n,  0  c  p.  ,  p„  cl  and 
3  j 


y  Pc  *  I  Pn  "I  and  where  we  are  assuming  that  S(w) 

3-1  3-1  "j 

and  N(w)  are  symmetric  about  the  origin  and  contin¬ 
uous  at  w-a, ,  3-1, 2,... m.  The  breakpoints,  a,,  and 

J  2  2^ 

and  the  fractional  powers,  p-  a.  and  p  o„,  are  as- 

S3  S  N 

sumed  known.  The  mth .breakpoint  Is  taken  to  be  In¬ 
finity  unless  the  bandwidth  of  the  signal  is  also 
These  classes  of  PSD's  are  termed  p-polnt 
spectral  classes,  and  are  useful  because  the  frac¬ 
tional  power  Is  an  easily  measured  parameter  of  the 
process . 

Intuitively,  It  seems  that  If  the  only  know¬ 
ledge  about  the  signal  and  noise  PSD's  involves  the 
fractional  powers,  then  the  least  favorable  signal 
and  noise  PSD's  can  be  anything  so  long  as  their 
ratio  Is  piecewise  constant.  The  previously  men¬ 
tioned  Image  processing  application  hints  at  this 
solution.  Also,  previous  work  in  robust  filtering 
[1,3]  has  shown  that  the  least  favorable  PSD's  tend 
Co  look  as  much  alike  as  possible.  These  Intuitions 
are  stated  more  rigorously  In  the  following  theorem. 


Theorm  1;  For  any  signal  and  noise  PSD's,  ,S(w) 
and  N(w),  which  are  members  of  the  general  p-point 
classes,  X  and )( ,  respectively,  the  robust  filter, 
Hjj(w) ,  satisfying  (18)  is  given  by 


Sj.(w) 

“r^”^  “  S^(w)+N^(w) 


(20) 


where  the  least  favorable  PSD's,  S^(w)  and  Wj^(v), 
are  such  that 


Sjj(w)  •  kjNjj(w),  wElaj_j^,aj)  (21) 


where 


P- 


S(w)dw  Pg  Og 

kj  -  3-1.2, ...  ,m 


r 


(22) 


N(w)dw 


"3 


•3-1 


Proof ;  The  right-hand  side  of  (18)  is  true  by  the 
definition  of  8|^(w)  as  the  optimum  filter  for  S^(w) 

and  N|^(w).  Thus,  It  remains  to  show  that 
e(S;N;Hg)  s  e(Sg,N^:H^)  for  all  (S,N)e1xK. 

The  appropriate  MSE's  are  calculated  as 

•<W'V-7  J, 

3-1  aj.i  3  3 

and 

a(8,N;V7  I  ?^(S(w)(^)^-HI(w)(^)^)dw 

-R  It  l+kj  l+kj 


But  since  S(w)  and  S|^(w)  both  belong  to  jt  and  K(w) 
and  Il^(v)  both  belong  toMi  both  quantities  reduce 
to  the  sane  result  and  we  obtain 

2  2 

e(S.N;H^).e(S^.N^;H^)  -  f  -1  ^  i  (jj) 

and  the  assertion  Is  proved.  Thus,  once  the  ro¬ 
bust  filter  Is  constructed  It  gives  the  same  level 
of  performance  for  all  signal  and  noise  PSD's  In 
the  general  p-polnt  class  for  which  It  was  defined. 
Contrast  this  with  the  performance  of  the  optimum 
Wiener  filter,  H^(w),  which  may  degrade  consider¬ 
ably  with  snail  deviations  from  the  assumed  PSD’s. 

It  la  instructive  to  consider  the  statement 
made  earlier  that  the  least  favorable  PSD's  tend 
to  be  "closest"  In  some  sense.  We  will  use  the 
concept  of  distance  measures  [7,8]  to  Justify  this 
Statement.  We  have  the  following  lemma. 

Ltamt  1;  The  least  favorable  PSD’s,  Sj^(w)  and 
H^<w)  ,  defined  as  In  Theorem  1,  minimize  the  dist¬ 
ance  measure,  d(S,N),over  the  p-polnt  spectral 
classes,  iS- and M ,  where  d(S,N)  is  defined  ast 

A  ? 

d(S.ll)  -  Ej,[ca)]  -  2/  H(w)ca)dw  (2A) 

where  C(*)ls  any  real,  continuous,  convex  function 
on  (o,«>)  and 

Proof;  Assuming  d(S,N)  Is  finite,  we  must  show 
that  the  least  favorable  pair  minimizes  the  dist¬ 
ance  measure,  l.e., 

dCS^.R^)  i  d(S,»)  for  all  (S,N)clxM  (25) 


How, 


l^re  we  are  assuming  that  0  <  Sj^(w)  ,N|^(w)  <  •>. 
But,  for  convex  functions,  C’’(w)  a  0,  ■<>  that 


S  S  ^R 

C(|)-C(J)RC’(5j)(|-p^) 

Equation  (26)  then  becomes  (after  some  manipula¬ 
tion) 

^d(S,H)-d(S^,H^)la  *  [/^(H-N^){C(jj^)-C’(^)(5^)}dw 

®  R  R 


S 

+  r  (S-Sj^){C'0)dvj  (27) 

*j-i  * 

However,  from  Theorem  1,  the  ratio  piece- 

wise  constant  so  that  (27)  becomes 

i{d(S,H)-d(S^.N^)]  > 

I  [Aj  J  (N-Ng)dw  +  B  J  (S-Sg)dw] 

‘j-1  *3-1 


where  A^  and  Bj  are  constants.  But  since  S(w)  and 
Sj^(w)  both  belong  to  and  N(w)  and  Hj^(w)  both  be¬ 
long  toN,  the  right-hand  side  is  zero  and  we  ob¬ 
tain 

d(S,N)  >  d(Sj,Nj^)  for  all  (S,N)eAxK 


the  desired  result. 


Example;  A  simple  (although  somewhat  contrived)  ex¬ 
ample  seirves  to  illustrate  the  performance  invari¬ 
ance  of  the  robust  filter.  Assume  the  nominal  sig¬ 
nal  and  noise  PSD’s  are  given  by 


S^(w) 


|w|  t  1 
otherwise 


l-|w| ,  |w|  s  1 
0,  otherwise 


2  2  1 

with  total  signal  and  noise  powers  Og'Ojj"  27* 

will  assume  that  these  PSD’s  are  really  "guesses" 
based  in  part  on  power  measurements  made  at  the  ar¬ 
bitrarily  chosen  p-point  locations  a^*l/2  and  a2~l 

(the  signal  bandwidth  value  which  has  been  assumed 
known).  The  optimum  filter  for  these  nominal  PSD’s 
is  obtained  as 


H„(w) 


c  1 


with  the  corresponding  minlaaim  HSE,  e(S^,H^;H^)  ■ 
0.0717.  The  robust  filter,  Hj^(w) ,  designed  using 

only  knowledge  of  the  fractional  signal  and  noise 
powers  Is  obtained  as 

H_(w)  -  ^  0  ^  f  1/2 

V’''  \  2/3,  1/2  S  lw|  <  1 

with  the  corresponding  MSE,  e(S^,N^;Hj^)-0.0743  (a 

degradation  of  only  '3.6Z).  Suppose,  now,  that  the 
actual  signal  PSD  (let  N(w)>N^(w))  Is 


S(w) 


{5: 


0  s  |w|  <  1/4,  1/2  i  |w|  <  3/4 
otherwise 


Ho(w). 


Using  the  previously  derived  optimum  filter, 

the  performance  degrades  '13Z  to  a  value  of 
e(S,Ho:Ho)*'0.0811  while  the  robust  filter  maintains 

the  ssme  performance  level,  e(S,N;H|^)>0.0743.  Ob¬ 
viously,  for  any  deviation  in  S^(w)  and  H^(w),  so 

long  as  the  resulting  PSD’s  are  In  the  general  clas¬ 
ses,  ^andM,  the  robust  filter  will  provide  the 


t  w  dependence  has  been  dropped  for  convenience. 


MM  level  of  perfonunce  ehlle  the  optlnum  filter 
pcrforeence  aay  degrade  severely. 

Up  to  this  point  we  have  not  confronted  the 
gnestton  of  how  to  choose  the  location  of  the  break¬ 
points.  Obviously,  if  nothing  but  the  fractional 
powers  are  known  it  Is  quite  arbitrary  where  we 
place  the  breakpoints.  However,  if  sone  additional 
knowledge  about  the  signal  and  noise  characteris¬ 
tics  Is  available,  we  can  use  this  knowledge  to 
better  choose  the  breakpoint  locations.  For  exam- 

If  "guesses"  are  available  for  Che  signal  and 
noise  PSD’s  we  can  use  the  OPC  filter  as  a  guide 
for  designing  the  robust  filter. 
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TABLE  1.  OPC  filter  parameters  and  corresponding 
HSE's  for  ButCerworth  PSD's  (MSB's 
using  Hq(w)  Is  given  for  comparison)  [9] 


p_ 

n 

*S 

m 

“j 

^1 

MSE 

MMSE(Ho(w)) 

3 

1 

1 

2 

1.392 

8:gi3 

0.4030 

0.3934 

3“ 

1 

1 

3 

1.184 

0.627 

0.3972 

0.3934 

1.681 

0,352 

0.039 

3“ 

1 

.5 

2 

0.670 

0.718 

0.063 

0.3085 

0.2939 

T" 

T 

1 

2 

1.103 

0.651 

0.022 

0.3557 

0.3499 

9" 

T" 

1 

3 

1.046 

0.655 

0.3531 

0.3499 

1.165 

0.351 

0.009 

9~ 

1 

1 

4 

1,038 

0.655 

0.3531 

0.3499 

1.153 

0,396 

1.810 

0.036 

0.001 

T 

1 

2 

'T~ 

2.277 

0.574 

0.022 

0.4291 

0.4084 

J~ 

-y- 

1 

2 

1.243 

0,529 

0.177 

0.4821 

0.4775 

T- 

“4* 

1 

2 

1.249 

0.471 

0.4821 

0.4775 

0.830 


